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CN ! Abstract 

We show that, in cosmological microlensing, corrections of order v/c ~ AA/A, 

to the deflection angle of light beams from a distant source are not negligible and 

r — . J that all microlensing quantities should be corrected up to this order independently 

CO ■ of the cosmological model used. 
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In the last decade, gravitational lensing has become one of the most powerful tools in 
astrophysics and cosmology in studies about the presence and mass distribution of dark 
matter in the Universe [Q , || , flU . It affords, in principle, estimates of the gravitational 
> ! mass of all large scale structures, from galaxies to super clusters, and, in the specific ap- 

plication called microlensing, it can be used to search for Massive Astrophysical Compact 
Halo Objects (MACHOs) ||. These objects are considered to be the main constituents 
of the dark halo of spiral galaxies (of our Galaxy, in particular) and, theoretically, could 
range in mass from 10~ 8 -j- 1O 2 M . MACHOs could therefore represent planets, brown 
dwarfs, or massive black holes ||. The fundamental problem is how lensing by a point- 
like mass can be detected. Unless the lens is very massive (M > 10 6 M Q ), the angular 
separation of two images, usually produced by a point-like lens, is too small to be resolved. 
The angular separations of images are, in fact, of the order ~ 10~ 3 -7- 10~ 6 arcsec, hence 
the term microlensing. However, even when detecting multiple images is impossible, the 
magnification can still be seen if the lens and the source move relative to each other. 
This motion gives rise to a lensing-induced time dependence of the source luminosity 0. 
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Since the effect was first observed for the quasars QSO 2237+0305 and QSO 0957+561 
0,11, one must now distinguish galactic microlensing from extragalactic or cosmological 
microlensing. In the first case, the light sources are stars and the angular separations in- 
volved are ~ 10~ 3 arcsec. In the second case, the sources are very distant quasars and the 
angular separations involved are ~ 10 _6 arcsec. In both instances the term "microlens- 
ing" is used. In this work we are exclusively concerned with cosmological microlensing. 
The principle of microlensing is quite simple. If the closest approach between a point 
mass lens and a source is equal to or less than Be, the Einstein angular radius, the peak 
magnification in the lensing-induced light curve results in a brightness enhancement (e.g. 
~ 0.3 magnitudes), which can be easily detected. As discussed below, the Einstein an- 
gular radius 9e is a feature of the lens-source system and provides the natural angular 
scale required to describe the lensing geometry. It gives, in fact, the typical angular sep- 
aration among the single images, while for axisymmetric lens-source-observer systems, 
it gives the aperture of the circular bright image, called Einstein ring. Geometrically, 
the Einstein ring, can be defined in any case, even when a luminous circular image is not 
produced. However, sources which are closer than 9e to the optical axis experience a 
strong lensing effect and are hardly magnified, while sources that are located well outside 
the Einstein ring are weakly magnified. Therefore, for a lot of lens models, the Einstein 
ring represents the boundary between the zones where sources are strongly magnified or 
multiply-imaged and those where they are softly magnified or singly-imaged ||. 

The first detection proposal || consisted in monitoring millions of stars in the Large 
Magellanic Cloud (LMC), or in the bulge of the Galaxy and in searching for the cor- 
responding magnifications. By detecting a sufficient number of events, one could then 
map the distribution of (dark) stellar-mass objects in the halo of the Galaxy (due to 
the fact that LMC is near us and the halo of our galaxy is in between) or between the 
Solar System and the bulge of the Galaxy. The two alternatives require some care in the 
selection of distances between source and observer. In fact, the distance between the Sun 
and the center of LMC is ~ 55Kpc while the distance between the Sun and the Galaxy 
bulge is ~ 8.5Kpc. This difference in size yields Einstein radii for the selected sources 
which could differ by about one order of magnitude. Furthermore, the Galaxy haloes are 
supposed to extend for ~ 50Kpc so that the zone where MACHOs can pass is very large. 
Both approaches can however be used for galactic microlensing and for r^ ~ 1 + 10AU, 
source-lens-observer distances D ~ 1 -f- 50Kpc and MACHO velocities v ~ 100 -f- 500Km 
s _1 microlensing effects may be observable @]. 

In the proposal of Ref. |||, distinguishing the intrinsic variable stars (which are very 
numerous in a normal galaxy) from the lensing-induced variables is potentially a serious 
problem. Fortunately, the light curves of lensed stars have distinctive features which can 
be used to separate induced variability from intrinsic variability. For instance, the light 
curves are symmetric in time and chromatic effects are absent because light deflection 
does not depend on the wavelength. On the contrary, intrinsic variables have asymmetric 
light curves. Furthermore, magnification produces in this case strong chromatic effects. 
The probability of seeing microlensing events depends on the optical depth, which is the 



probability that at any instant in time a given source be within the angle 9e of a lens. 
The optical depth is the integral over the number density n(D i) of lenses times the area 
enclosed by the Einstein ring of each lens, i. e. 
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where dV is the volume of an infinitesimal spherical shell of radius D \ which covers a 
solid angle Q. Eq.fll]) may take a very simple form if the sources are distant and compact. 
Then sources and lenses have angular sizes smaller than Be- 

From the cosmological point of view, microlensing plays a fundamental role in the 
determination of the density parameter Q Q and the cosmological constant A. It has been 
frequently argued that a significant fraction of the dark matter in the Universe may be 
in the form of compact masses which could induce lensing phenomena. For example, 
let us consider an Einstein-de Sitter Universe with a certain constant comoving number 
density of point lenses of mass M. Let Qm be the density parameter. The optical depth 
for lensing of sources at redshift z s can be shown to be 



t(z v ) = 3Q 
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which is t(z s ) ~ Vl m z 2 s JAl, for z s <C 1, and t(z s ) — 0.3f2 M , for z s = 2. Hence the number 
of lensing events in a given source sample measures the cosmological density of compact 
objects directly. These results depend strongly on the red-shift ||. In calculating the 
probability of lensing it is important to allow for various selection effects. Lenses magnify 
the observed flux, and therefore sources which are intrinsically too faint to be observed 
may be lifted above the detection threshold. At the same time, lensing increases the 
solid angle within which sources are observed so that their number density in the sky is 



reduced |TIJ . If the faint sources are numerous, the increase in source number due to the 
apparent brightening outweighs their spatial distribution, and the observed number of 
sources is increased by lensing. This magnification bias JTT|] can increase the probability 
of lensing for bright optical quasars. Besides, the strong dependence of the optical depth 
on the red-shift suggests that the simple weak field and slow motion approximations 
(normally used in lensing theory) be applied with care. In particular, we show below 
that discarding terms of the order c~ 3 in the lens equation is not justified in the high 
red-shift regime. In general, it is not always correct to consider only the term $/c 2 in 
the refractive index n because the propagation of light in a gravitational background 
depends on the mass distribution and, frequently, the ray approximation does not work 

0,0- 

In this letter we show explicitly that some effects ensue from including higher order 

approximation terms in the calculation of the light deflection by a standard lens (in 

particular a point-like lens). 

In general relativity the weak field approximation is defined by 

9nA x ) = Vtw + h^ix), \htj, v (x)\ < 1. (3) 



If the distances are small with respect to the Hubble distance c/Hq, one can neglect 
the curvature and the expansion of the Universe and the stress-energy tensor for perfect 
fluid matter is given by 

T" v = (p + pc 2 )u»u v - pg^ , (4) 

which, in the approximation \v\ <C c and p <C pc 2 , reduces to the components 

T 00 ~ pc 2 , T 0j ~ pcv j , T ij ~ pv l v j . (5) 

By using the Einstein equations, one finds 

T _ r2l 8ttG 21 87tG 2 167rG i , . 

V Zioo = — j-P. V ftij = — j-OiiP, V /ioi = —bjipv , (6) 

where V 2 is the usual Laplacian operator. In arriving at Eqs.(^|), use has been made of 
the harmonic condition 

<T r£, = o . (7) 

The integration of Eqs.ffl) yields 
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h 00 = 5- / r=; =77 a x , rt-y = 5-% / — —d x , ft j = -5-% / t^ ^t" £ . 
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(8) 
Up to leading order in u/c, the metric is determined by the gravitational potential 



$( x ) = -G ( —^—d 3 x' (9) 
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and by a potential V^ 

^ = -G/-^-rfV. (10) 
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From Eqs. (§) and (|8[)- (|T0"1) , one gets 

rfs 2 = (i + ^ c 2 rft 2 - ^-M^cdtdx 1 -(l-^\ S l ,dx i dx j . (11) 
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By calculating the affine connection related to the metric (|i~l|), one also obtains the 
geodesic equations 

x a + r^afx" = (12) 

where the dots indicate differentiation with respect to the affine parameter. These are 
the only ingredients necessary to derive the gravitational lens equation for a beam of light 
rays propagating in a weak gravitational field up to and including terms proportional to 

yl 

Equations ( |T2"| ) yield, to leading order in v/c, 
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Since for a light ray ds = da = 0, Eq. (JllD gives, to order c 

4V' , / 2$\ 
cdt = — — dx + I 1 -J rfZ eud , (15) 

where c?Z 2 ud = 5ijdx l dxi is the Euclidean length interval. Squaring ( |i"5|) and keeping 
terms to order 1/c 3 , one finds 

/ 4-<I>\ 81/' 
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Inserting ( jig ) into Q14]) , one gets 
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From ( |T7| ) we derive the equation of the gravitational lens. Let us consider l euc i as a 

parameter. In the approximations used, da ~ rfZ euc z and Eq. fli"3|) can be written in the 

form 

d fcdt 2<&cdt A ,dx l \ 

da \ da c 2 da c 3 da J 

from which 

cdt / 2$\ 4,cfe z 
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The affine parameter can be chosen to make the constant in Eq. (|1^) unity. Eqs. (|15|) and 
(|19|) imply that, to lower order in v/c, 

dleud = 1 

Therefore, for weak gravitational fields, the vector tangent to the trajectory of a light 
ray can be expressed as 

dx k dx k , , 

(21) 
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It is then possible to introduce the vector 



dx k 
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and Eq. (|17|) can be recast in the form 
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using the relation 
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in Eq. (g3j) one finds 
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The first term is the component of the gradient of <l> orthogonal to the vector e, i.e. 
Vj_ = V — e(e ■ V), while the second term, is the gravitomagnetic-term. Eq . (|25"D reads 
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where e is a unit vector. The deflection angle a of a light ray propagating in a weak 
gravitational field is given by a = e{ n — e ou t- In the general case, one finds from Eq.( 
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where the first term behaves as $/c 2 and the second as $t> /c 3 . 

Alternatively, the equation of a lens can be recovered by using Fermat's principle. In 
fact, in classical optics, the light rays follow trajectories which minimize the optical path 
Jndl. One can solve the equation ds 2 = with respect to the temporal coordinate dx°. 
The trajectory of a light ray is the extremal of the integral 
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with dleuci — Sijdx l dxK The integrand in Eq. (|2~8|) is the refractive index n of a light ray 
propagating in a gravitational field. In the weak field approximation, one can use Eq. 
(0) to get 

5 I f ^-dx l +(l-^L) JSijdxtdxi 1 = 0, (29) 
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which gives the trajectories of light rays. It therefore follows that 
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By using the Euler-Lagrange equations and using the previous results, one easly obtains 
the deflection angle a given by Eq . (|27|) . For small deflection angles and weak gravitational 
fields, which are the regimes of practical interest, the true position of a light source on 
the sky relative to the position of its image (s) can be defined. The lens equation can be 
recast in the form 

9- 0s = - (^r-) u0) = & , (31) 



D. 



OS 



where 9 is the position(s) of image(s) with respect to the optical axis, 9 S the position 
of the source and a is the displacement angle. _D/ S and D os represent respectively the 
distance between the lens and the source and that between the observer and the source. 
In general, a given image position always corresponds to a specific source position, 
whereas a given source position may correspond to several distinct image positions. For 
point-mass lenses, the geometry of the system is simplified and we need not use the full 
vector equation (pH|). In this case, the deflection angle to order c~ 3 is given by 

A AGM 8GMv 

a ~ — , (32) 

ct c 6 r 

where tq is the impact parameter. From the definition of red-shift, AA/A = v/c, it 
follows that the deflection angle can be rewritten as 

^ AGM 8GMAX AGM ( 2AA\ , , 

If AA/A is not negligible, the microlensing quantities can exhibit an interesting behaviour. 
By writing r = 9D i, the lens equation for a point-mass lens takes the form 



which can be rewritten as 
where 



9 2 - 9 S 9 - 6% = , (35) 

is the square of the Einstein angle defined by r E = 9 E D oi . 9e depends on the distances 
involved and the mass of the deflector. The symbol M(< r E ) signifies that the mass of 
the lens must be be contained inside a sphere of radius r E . 

Before solving the algebraic Eq.(^), the important parameter magnification must be 
discussed. Gravitational lensing preserves the surface brightness of a source. Then the 
ratio of the solid angle dfli covered by the lensed image to that of the unlensed source dQ s 



gives the flux amplification (magnification) due to lensing. This is given by the Jacobian 
of the transformation matrix between the source and the image(s) 
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(37) 



If there is more than a single image, the total magnification is the sum of all image 
magnifications. Considering, as we do, a gravitational point-mass lens which is axially 



symmetric with respect to the line-of-sight, we can use the scalar angle (|32| ) and apply 
Gauss's law for the total flux. The light deflection reduces to a one-dimensional problem 
and Eq.([37]) becomes 

_ OjdOj 

6 s d9 s 

which can be easily applied 0. Let us now solve Eq . (p5|) . We get 



o±-^±\l^ + o 2 E ; 
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from which we see that 



e, = o 



±e t 



(39) 



(40) 



Because of ([36]), the position of the images is therefore shifted by AA/A ~ v/c. Eqs. 



(fl0|) indicate that we must expect at least two images from the same source. These lie on 
the same plane of the source. It also follows that all quantities where 9e and 9± appear 
are modified by AA/A. As discussed above, it is in general difficult to separate the two 
images and the result is a luminosity enhancement of the source. The magnification 
corresponding to Eq.(|39|) is 

which shows that when 9 S is zero, the magnification becomes singular. Physically, this 
means that when the optical system source-lens-observer is aligned, we can get a huge 
magnification. The total amplification due to both images is 
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where x = @s/9e- It immediately follows that 9 S < 9e — ► [i> > 1.34, which is the 
condition for the magnification inside the Einstein ring: a magnification \i ~ 1.34 corre- 
sponds to a magnitude enhancement Am ~ 0.32 as required in microlensing experiments. 
In other words, when the true position of a light source lies inside the Einstein ring, the 
total magnification of the two images amounts to /i > 1.34. This means that the angular 
cross section necessary in order to have significant lensing (i.e. fi ~ 1.34 and Am ~ 0.32), 
is equal to tt9 e which, from (1361), is proportional to the mass M of the deflector and to 
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the ratio of the distances involved. We can now calculate the optical depth in the case 
of randomly distributed point-mass lenses. It is possible to estimate the frequency of 
significant gravitational lensing events from the observations of distant compact sources. 
This is equivalent to considering optical systems with angular sizes much smaller than 
Oe ■ In this situation, the magnification of a compact source is equal or greater than 
1.34 (since 9 S < 9e) and the probability P of significant lensing for a randomly located 
compact source at a distance D os is given by 

„ tt91 ( D u \ (GM\ / AA' 



(m 1-2^ , (43) 
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where use has been made of flSTp. Eq. (H3) is linear in M and therefore holds true 



also when several point-mass lenses are present because the masses can be summed up. 
Assuming a constant density for the lens(es) and a static background (this last assumption 
surely holds for galactic distances), and averaging over the distances Di s ,D i, D os , the 
probability (|43|) can be interpreted as the optical depth r for lensing [jl4|| , flBH , P| . 



The considerations given above show that the c -3 corrections affect all microlensing 
quantities, i.e. a, ##,//, r. Some estimates are now in order. Let us consider an ob- 
ject receding from us with a velocity v ~ 10 3 Km/s. SinceAA/A = v/c, the correction 
is one percent and could be measured with some difficulty If the recession velocity is 
v ~ 10 4 Km/s, the correction is ten percent and could in all likelyhood be observed. Far 
objects like quasars can easily possess similar recession velocities. It is therefore conceiv- 
able that by carrying out observations with precision higher than presently available it 
would be possible to observe the effects discussed in the case of cosmological microlens- 
ing. Essentially, they should consist in a shift of image positions and in a reduction of 
the amplification curve. Extremely precise observations could, in our opinion, detect the 
effects. In a forthcoming paper, we will apply these results to other lens models such as 
the isothermal sphere or the disk of galaxies. 
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